Holon-pair boson theory based on the U(l) and SU(2) slave-boson approaches to the 

t-J Hamiltonian 
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To supplement our recent brief report on the theory of holon-pair boson approach to the t-J 
Hamihonian [S.-S. Lee and Sung-Ho Suck Salk, Phys. Rev. B 64, 052501(2001)], in this paper we 
present a full exposure to the theory, detailed physical implications and predicted various physical 
properties of high Tc cuprates. We discuss the significance of coupling (interplay) between the 
spin and charge degrees of freedom in the Heisenberg interaction term of the t-J Hamiltonian. We 
discuss its importance in causing the arch-shaped superconducting transition temperature Tc and 
the pseudogap (spin gap) temperature T* tangential to Tc in the overdoped region in the observed 
phase diagram of high Tc cuprates. A universal parabolic scaling behavior of T* /Tc (or Tc/T*) 
with hole doping concentration is predicted in agreement with observations, indicating that there 
exists correlation between the pseudogap (spin gap) phase and the superconducting phase through 
antiferromagnetic fluctuations. Our proposed holon-pair boson theory is shown to be self-consistent 
in that it not only yields the arch (dome) shape structure of Tc but also reproduces numerous other 
physical properties such as superfluid weight, bose condensation energy, spectral function, optical 
conductivity and spin susceptibility, including their temperature and doping dependence. 



I. INTRODUCTION 

Since the discovery of high Tc cuprates [1], the pseu- 
dogap phase and /superconducting phase have been re- 
vealed in various measurements [2] such as the angle re- 
solved photoemission spectroscopy (ARPES) [3-6], tun- 
neling spectroscopy [7-9], nuclear magnetic resonance 
(NMR) [10-13], magnetic susceptibility [8], inelastic neu- 
tron scattering [14], optical conductivity [15-17], Ra- 
man scattering [18,19], resistivity [20] and specific heat 
[21-24]. The arch shaped bose condensation (super- 
conducting transition) temperature and the tangen- 
tial appearance of the pseudogap temperature T* to 
in the overdoped region are commonly observed in the 
phase diagram of both high Tc cuprates La2-xSrxCu04 
and Bi2Sr2CaCu208+5 [8,3,10-13,18,19,21]. Most re- 
cently, from the Nernst effect measurements in the high 
Tc cuprates, Bi2Sr2CaCu20s-i-s Wang et al. showed that 
in agreement with the observed trend of the ARPES gap 
[5] the Cooper pairing potential monotonically decreases 
with increasing hole concentration [25]. The pseudogap 
also decreases with increasing hole concentration. Paying 
attention to the ratio of the observed spin gap (pseudo- 
gap) temperature T* to the superconducting tempera- 
ture Tc, Nakano et al. [8] observed that there is a good 
correlation between T* and Tc, and that there exists a 
universal scaling behavior with hole doping, that is, the 
sample independence in the ratio T* /Tc (or Tc/T*) for 
the entire range of hole doping concentration for both 
La2-xSrxCuOi and Bi2Sr2CaCu20s+s- This hints that 
superconductivity is affected by (correlated with) spin 
fluctuations or spin singlet pair excitations. The generic 
features of peak-dip-hump structures have been observed 
in ARPES (angle resolved photoemission spectroscopy) 
[5,6], STM (scanning tunneling measurement) [7-9] and 



optical conductivity measurements [15], suggesting a pos- 
sibility of common origin. 

From the very beginning of high Tc superconductiv- 
ity, Anderson [26] has stressed that the t-J Hamiltonian 
is an appropriate model for the study of high Tc super- 
conductivity. Resonating-valence-bond (RVB) mean field 
approach to the t-J Hamiltonian [27] was introduced to 
unsatisfactorily predict a trend of monotonically decreas- 
ing superconducting transition temperature with increas- 
ing hole doping. Later by applying the U(l) slave-boson 
theory of the t-J Hamiltonian [28-32] , decreasing pseu- 
dogap temperature with hole doping concentration was 
predicted and identified to be the spin gap temperature. 
Superconductivity was proposed to arise as a result of 
bose condensation of composite particles made of the 
spinon singlet pair and single holons. Later, Wen and 
Lee [33] proposed an extended SU(2) slave-boson the- 
ory to impose symmetry not only for the spinon but for 
the holon for hole doped systems. In theories which pay 
attention to the bose condensation of single holons, the 
predicted bose condensation temperature Tc showed a 
linear increase with hole concentration x. In these the- 
ories Tc is scaled by the kinetic energy of doped holes 
(of the order of hopping energy tx) and the pseudogap 
temperature T* by the Heisenberg coupling J. Accord- 
ingly Tc was shown to increase linearly with increasing 
X. These theories fails to predict correlation between Tc 
and T* and thus the experimentally observed universal 
behavior of T* /Tc can not be explained. Gimm et al. 
[34] proposed a theory of holon-pair bose condensation 



term in the Heisen- 

i„ 



by paying attention to the 

berg interaction term of J{Si ■ Sj — ^riirij). As a result, 
the bose condensation (superconducting transition) tem- 
perature scales with the Heisenberg coupling J but failed 
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to predict the arch (dome) shaped superconducting tran- 
sition temperature in the observed phase diagram. 



To remedy this failure, lately we proposed an improved 
version of holon-pair boson theory [35] which takes into 
account coupling between the charge (holon) and spin 
(spinon) degrees of freedom in the slave-boson repre- 
sentation of the Heisenberg interaction term in the t-J 
Hamiltonian. Theories which neglect the charge con- 
tribution in treating the interaction term are expected 
to fail in providing superconductivity. This is because 
the Cooper pairs represent pairing of charged particles 
(— 2e or +2e). In such incomplete theories only the spin 
(spinon) paring is possible. On the other hand, in our 
theory both the charge and spin pairing are possible since 
the Heisenberg interaction term is represented by cou- 
pling between the charge and spin contributions. Such 
coupling is proven to cause the arch-shaped supercon- 
ducting transition temperature [35]. Further, correlation 
between the spin gap temperature T* and the supercon- 
ducting transition temperature Tc is predicted in agree- 
ment with observation [8] . The coupling between the spin 
and charge degrees of freedom results in salient predic- 
tions of various physical properties; they are the arch- 
shaped superconducting temperature [35], bose conden- 
sation energy [36], spectral function [37], the peak-dip- 
hump structure in the optical conductivity [38], the uni- 
versal scaling behavior of T* jT^ with hole density x [37], 
the superfluid weight vs. Tc and the doping and temper- 
ature dependence of superfluid weight [39] . 



It is reminded that our approach is pivotally different 
from other proposed U(l) and SU(2) slave-boson theories 
[28-31,33] in that the couphng of the spin and charge de- 
grees of freedom naturally appears in the slave-boson rep- 
resentation of the Heisenberg interaction term and that 
physical spin-charge separation does not appear even in 
the mean field Hamiltonian. As a result of the coupling, 
this theory predicts the spin-gap (pseudogap) tempera- 
ture T* tangential to in the overdoped region and a 
universal behavior of the ratio T* /Tc with x. Realiz- 
ing that the pseudogap is identified as the spin gap from 
NMR measurements [10-13] and that Tc is observed to 
be correlated with T* [8], we discuss how the spin sin- 
glet pair excitations (spin gap) affect the high Tc super- 
conductivity. Our earlier presentation of this holon-pair 
boson theory [35] is excessively brief. In this paper a full 
disclosure to the details of the holon-pair boson theory 
will be made with emphasis of physical implications, and 
theoretical predictions, highlighted. By means of self- 
containment we show the self-consistency of our theory 
by revealing that not only the phase diagram but also 
various other physical properties are well reproduced in 
agreement with observations. 



II. INTERPRETATION OF THE HEISENBERG 
INTERACTION TERM IN THE U(l) 
SLAVE-BOSON REPRESENTATION 

The t-J Hamiltonian of interest is given by, 

H = -tY^ {c\^Cja + c.c.) + J ^{Si- Sj - ^rurij) 

<i,j> <i,j> 

-/^X!4<7Cicr (1) 
i 

and the Heisenberg interaction term is rewritten 
Hj = J ^ (Si • Sj - ^riiUj) 

<i,j> 

= - ^ XI ^ 4]4l')('^3T^^l " CjiQi). (2) 

<id> 

Here t is the hopping energy and Sj, the electron spin 
operator at site i, Si = \c\c^(t afjCijs with cTq/j, the Pauli 
spin matrix element, rij is the electron number operator 
at site i, Tii = cl^Cia- n is the chemical potential. 

In the U(l) slave-boson representation [28-31], with 
single occupancy constraint at site i the electron annihi- 
lation operator Cia is taken as a composite operator of the 
spinon (neutrally charged fermion) annihilation operator 
fia and the holon (positively charged boson) creation op- 
erator h\, and thus, Ci^ = fiab\- Rigorously speaking, it 
should be notc^d that the expression Ci„ = h\ji„ is not 
precise since these operators belong to different Hilbcrt 
spaces and thus the equality sign here should be taken 
only as a symbol for mapping. Using Cj^ = !i<j^i and 
introducing the Lagrange multiplier term (the last term 
in Eq.(3)) to enforce single occupancy constraint, the t-J 
Hamiltonian is rewritten, 

H = -tYl {{flh){b]fj.) + c.c) + Hj 

<i,j> 

i 

-iJ2\{blh + flfi^-l) (3) 

i 

with the Heisenberg interaction term, 

<i,j> 

(4) 

The first term represents hopping of a spinon from site 
j to site i and of a holon (positively charged boson) 
from site i to site j. In the slave-boson representation 
a charged fermion (electron or hole) is taken as a com- 
posite particle of a 'spinon' and a 'holon'. They can con- 
veniently serve as book-keeping labels to discern physical 
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properties or objects involved with the charge or spin de- 
gree of freedom (e.g., spin gap phase, spin singlet pairs, 
hole pairs, ...). With the single occupancy constraint, 
electron is allowed to hop from a singly occupied copper 
site i only to a vacant copper site j. A site of single 
occupancy in the Cu02 plane of high Tc cuprates phys- 
ically represents an electrically neutral site (net charge 
0) with an electron of spin 1/2 and the vacant site, a 
site of positive charge +e with net spin 0. In the slave- 
boson representation, hopping of an electron (a compos- 
ite of spinon and holon) from a singly occupied copper 
site (neutral site) j to an empty site (positively charged 
site with +e) i results in the annihilation of a spinon (a 
fermion of charge and spin 1/2) and the creation of a 
positively charged holon (a boson of charge +e and spin 
0) at site j while at the copper site i a composite of a 
spinon (fermion of charge and spin 1/2) and a nega- 
tively charged holon is created. It is of note that as a 
result of electron hopping the newly occupied copper site 
i in the Cu02 plane can, also, be labeled as 'spinon' since 
this is an electrically neutral (charge 0) site with an elec- 
tron of spin 1/2 and the vacant site j, 'holon' since this 
is a positively charged site with a single charge +e and 
the net spin as mentioned above. Thus in practical 
sense, there is no distinction between the two different 
cases above. At times, we will call the singly occupied 
site as 'spinon' and the vacant (empty) site as 'holon' as 
long as there is no confusion. This is because any site oc- 
cupied by a spinon is identified as an electrically neutral 
site occupied by a single electron with spin 1/2 and the 
site with a positive holon is a positively charged vacant 
site with spin 0. Thus physical spin-charge separation is 
not allowed. 

The Heisenberg interaction term, Eq.(4) shows cou- 
pling between the charge and spin degrees of freedom. 
Physics involved with the charge degree of freedom is 
manifested by the four holon (boson) operator bibjbjbl 
in the Heisenberg interaction term. Judging from the in- 
tersite charge coupling term jUirij present in the Heisen- 
berg interaction term Hj = J '}Z<i j>{^i ' ~ jUirij), 
it is obvious that this charge contribution can not be ne- 
glected in its slave-boson representation. It is to be noted 
that the Hubbard Hamiltonian contains repulsive interac- 



tion U between charged particles and is mapped into the 
t-J Hamiltonian Ht-j in the large U limit. The Coulomb 
repulsion, Ujii^rni = ^{riif + riii)'^ - ■j('^»T ^ "-u)^ ob- 
viously manifests the presence of both the charge (the 
first term) and spin (the second term) degrees of free- 
dom. Thus, under mapping the charge part of contri- 
bution naturally appears in the Heisenberg interaction 
term. 

Let us now take another look at the importance of the 
charge contribution. In general, uncertainty principle be- 
tween the number density (amplitude ) and the phase of 
a boson order parameter applies. As an example, ar- 
bitrarily large fluctuations of the number density fix the 
phase, or arbitrarily large phase fluctuations fix the num- 
ber density of the boson. The conventional BCS super- 
conductors of long coherence length meet the former clas- 
sification, and thus the phase fluctuations of the Cooper 
pair order parameter are minimal. For charged bosons, 
e.g., the Cooper pairs, the number density fluctuations 
refer to charge density fluctuations. For short coherence 
length superconductors such as the high Tc cuprate sys- 
tems of present interest, local charge density fluctuations 
exist and cause large phase fluctuations. Thus, both the 
charge and phase fluctuations need to be taken into ac- 
count to fully exploit the quantum fluctuations . 

Let us now consider the importance of the charge and 
spin fluctuations. In generally, coupling between physical 
quantities A and B is decomposed into terms involving 
fluctuations of A, i.e., {A— < A >) and B, i.e., {B— < 
B >), separately uncorrelated mean field contribution of 
< A > and < B > and correlation between fluctuations 
of A and B, that is, {A- < A >) and {B- < B >) 
respectively; AB = {A- < A >) < B > +{B- < B > 
)<A> + <A><B> +{A- < A >){B- < B >). 
Setting A — bibjbjjbl for charge (holon) contribution and 

B = (//|/j| - fyji){f]hkl - fjl-k) for spin (spinon) 
contribution, the Heisenberg coupling term, Eq(4) can be 
decomposed into terms involving coupling between the 
charge and spin fiuctuations separately, the mean field 
contributions and coupling (correlation) between fluctu- 
ations (charge and spin fiuctuations). Using such decom- 
position of the Heisenberg interaction term for Eq.(3), 
we write the partition function, 



Z = j DbDf^DfiDXe-^^^- 



(5) 



where 



S[b, f, X] = dr bldrh + ^ fldrfi. + 



(7(1) 
t-J 



(6) 



with (3 = the inverse temperature and H^S]j , the U(l) symmetry preserved Hamiltonian, 



H^}'^ = -tY,{flf,J^b, + c.c.) 



<i,j> 



3 



E [((/lt/jT-4/]i)(/.T/a-/.i/n))Mi6t6t 

<i,j> 

+{b'^b,b]biy,flf^^^ flfj,){fnfii - fM 

+{bibjb]bl - {b,bjb]bl)) [ifyj^ - fyiM^M - fjiM - {{flifh - fhfhmfil - fnM) 

-f^T^flfiA^ + bib,) -iJ2x,{flf,^ + b% - 1). 



(7) 



III. U(l) MEAN FIELD HAMILTONIAN 

Noting that 6]] = % for boson, the intersite charge 
(holon) interaction term (the second term) in Eq.(7) is 
rewritten, 



fc=0 



J 



^{ifhfh - fyjoifnh - fM)bib,b]b\ 



(12) 

(13) 
(14) 



Here cr" is the unit matrix and a^''^'^ , the Pauli spin ma- 



-- < I Af. I >2 (l + b\bi + b]bj + blbtbjbi) , (8) trices. 

Combining Eq.(8) and Eq.(lO), we have 



with Ajj = fj\fii — fjifii, the spinon pairing field. The 
third term in Eq.(7) represents the intersite spin (spinon) 
interaction and is rewritten, 



J 



-<\A{A'>(l + blb, + b% + blb\b,b.. 



-2v^b,b]b])ifif], fy],mfii - fM 



(9) 



where Jp = J(l+ < bjb^ > + < b]bj > + < b\bib]bj >) 
or Jp = J(l — x)^ with a;, the uniform hole doping con- 
centration [40]. The fourth term in Eq.(7) is written. 



J 



^{bibjbibi){{fy]^ - fyj^mM - fM] 

J( + + + + \ ^J = -i;Y \^l?b\b%bi 

= ^(}^<b\b,>^<b\bi> + <b\b,b\bi>)<\A{.\^>. 2 Z-^i 

^^^^-^v E [^(4/jT-4/]i)(/^T/.l-/.l/iT) 
The intersite spinon interaction term in Eq.(9) is de- <«>J> 
composed into the direct, exchange and pairing channels \}_((t\ t \(t\ t 



ft ft 



+^ < \A{.f > (l+ < blbi > + < blbj > + < blbib]bj >) 

= < \A{.\' > blb]bjbi + I < \A{.\^ X blhb]bj > 

~ < |Af,f > [{bib,- < bib, >) + {bib,- < bib, >)] . (15) 

Collecting the decomposed terms Eq.(8) through 
Eq.(lO) in association with Eqs.(ll) through Eq.(15), we 
write 



J 



[31] 



fc=0 



^[E(/l^a/5/./3)(4^7^/.-^) - iflAfMlAsfjs)] +^ E < bib. X bib, > 



fe=l 

= Vd + Ve + Vp 



(11) _ J 



- 2 E [{b\bi- < b\h >) + {bib,- < bib, >)] , (16) 
with (j" = /, the identity matrix and fj^'^'"^, the Pauli <«J> 
spin matrices, where vd, ve and vp are the spinon inter- 
action terms of the direct, exchange and pairing channels where we considered < | Af^- [^ >= | Af^- [^ and ignored the 
respectively [35], fifth term in Eq.(7). 
I 



Hubbard Stratonovich transformation for the holon pairing term (the second term of Eq.(16)) leads to 



4 



and the saddle point approximation yields, 



Hi 



9 I y 



Ao;'|2-A-*(6,6,)-A-(6]6l) 



. Ob* 



(17) 



(18) 



where A^j* =< hihj > is the saddle point for the holon pairing order parameter A^^-. Since confusion is not likely to 
occur, we will use the notation A-'^ in place of A^j* for the saddle point. As are shown in Eqs.(12) through (14) the 

spinon interaction term is decomposed into the direct, exchange and pairing channels respectively. Proper Hubbard- 
Stratonovich transformations corresponding to these channels and saddle point approximation leads to the effective 
Hamiltonian, 



<i,j> 



<i,j> 



<»>j> 



+Y E EM)'-^i(/^'/^)] + YE(4/i.) 



<i,j> 1=0 



+4- E iblbj){b]h) 



<hj> 



4 y: \A{^^<bih><bib,> 



<i,3> 



- ^ \A{^f [(bib,- < bib, >) + (b]b,- < bib, >) 

<i,j> 

E 4/-(l + blh) - i ^ Xiiflfi. + blbi - 1), 



(19) 



where Af^- =< bibj >, Xij =< flfj^ + f/,bj >, 

=< fj^fn - fnh > and p'l =< IfU^fi > are 
proper saddle points. 

We note that p[ = \ < fja^f, >=< for 

/ = 1, 2, 3, /9° = i < 4/i<T >= i for / = for the contri- 
bution of the direct spinon interaction term (the fourth 
term). The expression {bljbi)(blhj) in the fifth term of 
Eq.(19) represents the exchange interaction channel. The 
exchange channel will be ignored owing to a large cost in 
energy, U ^ ^ [31] [33]. The resulting effective Hamil- 
tonian is 



(20) 



<i,j> 

4 El 



2\ \ b |2 



A/JV 



<i,3> 



ALI^ + -|x.,l 



1 

4J 



H is the holon Hamiltonian, 



<^,3> 

J 



- Af;(6,6,-)+c.c. 



<»j> 



-YiA{b\h-x), 



(21) 



(22) 



where iJ^'^ represents the the saddle point energy in- 
volved with the spinon pairing order parameter A-'', the 

holon pairing order parameter A'' and the hopping order where ix\ = + ^ 2j=i±s i±y lA^ P and , the spinon 



parameter x, 



Hamiltonian, 



5 



<i,j> 



<i,j> 



(23) 



where = — a;) + zAj. 

As can be seen from Eqs.(21) through (23), Eq.(20) 
reveals the importance of couphng between the spin and 
charge degrees of freedom, that is, couphng between the 
spinon pairing and holon pairing. Thus no spin-charge 
separation appears in the mean-field Hamiltonian above 
contrary to other mean field theories [28-31,33] which 
pay attention to the single-holon bose condensation. As 
can be seen from the second term in Eq.(22) which repre- 
sents holon pairing contribution it is expected that, owing 
to the coupling effect, bose condensation (or supercon- 
ducting phase transition) will occur only in the presence 
of the non- vanishing spin singlet pairing order, A-^ owing 
to the coupling effects mentioned above. Indeed, in high 
Tc cuprates superconductivity is not observed above the 
pseudogap (spin gap) temperatures T* where the spin 
singlet pairing order disappears. 



IV. U(l) FREE ENERGY 

The diagonalized Hamiltonian for Eq.(20) above is ob- 
tained to be (see Appendix A for detailed derivations), 



A 2 A 2 I A 2 „2 



fe,s=±l 



fe,s=±l 



where E^^ is the quasispinon excitation energy. 



(24) 



(25) 



with the spinon pairing energy (gap), Aq = Jp^k{T-^)Af, 
and E^^ is the quasiholon excitation energy. 



EL = W> 



(26) 



where the holon pairing energy, Aq = JAj^fe(r'')Ab and 
with (j) = 6, or r*, 



V7fcCOs2, 



i -I- ifl sin^ ( 



^ks 



2isx6(^), 



(27) 

(28) 
(29) 



with7/c = (cos fca;-|-cos fcj/) and (fik = {cos kx — cos ky) . X] 
denotes the summation over momentum k in the half re- 
duced Brillouin zone, and ,s = +1 and —1 represent the 
upper and lower energy bands of quasiparticles respec- 
tively. Here aksticxlsf) ^^'^ '^ksii 

'^ksi ') annihila- 
tion(creation) operators of spinon quasiparticles of spin 
up and spin down respectively, and PksiPlg), the annihi- 
lation(creation) operators of holon quasiparticles of spin 
0. e{g and e^^ are the kinetic energies for spinous and 
holons respectively. The minus sign (— A^) in the expres- 
sion of the holon quasiparticle energy ^(e — /z)^ — A^ 
arises as a consequence of the Bose Einstein statistics 
[41]. From the diagonalized Hamiltonian Eq.(24), we cal- 
culate the total free energy. 
Rewriting Eq.(24) as 



fe,s=±l 



^uw = E [4si4s^<^ks^ - ^ksialO + ElsPl^ks 

(30) 



with 



= NJA) [aI + x^) + NJj>(a) + y + J) 
+ E -Nxti^ + Nxn\ (31) 



fe,s=±l 

the partition function is derived to be, 

Z = exp(-/3if,) n (2cosh^)2(l-e-^^-)-i. (32) 

/s,s=±l 

Using the above expression, the total free energy is given 
by 

Fu(i) = NJA} (Ag + x^) + NJp [a] + Y + 1 
-2kBT ln{cosh{pElj2))- Nxn^ -2NkBTln2 

k,s=±l 



+kBT J2 ln(l-e-^<)+ ^ 



EL+f^' 



k,s=±l 



k,s=±l 



+ Nxff. 



(33) 



The set of uniform phase {9 = 0) for the hopping order 
parameter, d-wave symmetry [t^ = 7r/2) for the spinon 
pairing order parameter and s-wave symmetry (r** = 0) 
for the holon pairing order parameter is found to yield 
a stable saddle point energy for both the underdoped 
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and overdoped regions. There is another set of order pa- 
rameters which yield the same energy as the above one; 
27r-fliix phase {6 = 7r/2) for the hopping order param- 
eter, s-wave symmetry (r-'' = 0) for the spinon pairing 
order parameter and d-wave symmetry (r'' = tt/2) for 
the holon pairing order parameter. In both cases, the d- 
wave symmetry of the electron or hole (not holon) pairs 
occurs as a composite of the d-wave (s-wave) symmetry 
of spinon pairs and s-wave (d-wave) symmetry of holon 
pairs. Only at very low doping near half filling, the flux 
phase [31] becomes more stable. Thus, the phase of the 
order parameters of present interest are 6 = 0,Tf = tt/2 



and 



0. Then the d-wave symmetry of the electron 



or hole (not holon) pairs is a composite of the d-wave 
symmetry of spinon pairs and s-wave symmetry of holon 
pairs. Minimizing the free energy with respect to the am- 
plitudes of the order parameters %, A{, and Aj^, we obtain 
the self-consistent equations for the order parameters, 



dx 



dFu{i) 



dAf 




= 0, (34) 



ks 

INJAlA, 



ks 



+ ^ 



e/3sL -1 2j \ dAb 



0, (35) 



2iVJpA/ -^27VA/(A^ 



I tanh 



ks 



dAf 



ks 



1 2 



dA 



0. (36) 



For fixed niunb(>rs of spinon and holon at a p;i\'en hole; 



concentration, we obtain, for the chemical potentials, 
and yu'', 



dFu(i) 



I tanh 



PEL 



fe,s=±l 




Nx = Q, (37) 



fe,s=±l 



-I- 



-ks 



EL 



EL 



Nx = 0. 



(38) 



Using the five self-consistent equations of Eqs.(34) 
through (38), we determine x, Ab, A/, fj,-^ and fj,'' at 
each doping and temperature. Both the pseudogap tem- 
perature T* and the superconducting transition (bose 
condensation) temperature are determined to be the 
temperatures at which the spin gap Aq and the holon 
pairing energy (gap) Aq respectively begin to open. 



V. SU(2) ACTION FROM THE U(l) ACTION 

The t-J Hamiltonian is manifestly invariant under 
the local SU(2) transformation gi = e*"^'^ for both the 

fil\ _ ^ I fi] 



spinon and holon spinors with ( ^ 

Jil 



i2 



-fk 



m 



hi 

-A 



and 



bii 



9i 



hi 



satis- 



fyingc,T=6l/iT = (^'l>0)( f'/ 

hi 



(4,4) ( I 



and C; 



'ii 



blhi 



fii 

ibio)i !)\ 

J if 



fi: 



[33]. We intro- 



duce additional Lagrange multiplier terms involved with 
the constraints fj-^f}^ = and fiifi-[ = to write 



^■^Hflh-^+bibi-i)-iY,xjyi-iY,x;fiih 



it 



(39) 



in order to allow for SU(2) symmetry. Thus writing spinors V'ii = { j and = ^ j fo'' spinon, h'^ — 



for holon and the three-component Lagrangian multiplier field with a, 
the U(l) action in Eq.(6) can be rewritten as 



0(1) _ i\'+ix" 0(2) 
— 2 ' "j 



-\+X 0(3) 
2 ' "i 



iXi. 



-o E hih%h%hZ{fy],-fyimfii-fnfn) 



<ij> 



Y^flf^Mh + 1) - E^° • (^^^'-C + hV<ri^) 



(40) 
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Here the fourth term is the Heiscnbcrg interaction term, Hj — —-^ ^<_i,j> 
We rewrite the spinon part of the Heisenberg interaction, 



~'2^filfh ~ fitfjOifilfil ~ fjlfn) 



J 

fe=l 



(41) 



where = j and {fLa'^^fip) = \tr (*°^*?(£7'=)'^) [42]. Here denotes the transpose of the Pauh 

matrices for fc = 1, 2, 3. 



Realizing hi = gih^ and = 
interaction term is given by 



i2 



ft _ft 
Ji2 Jil 



— di [ j-t ft ), and using Eq. (41), the SU(2) symmetric Heisenberg 



fll -I. 



it 



<hj> 
J 



- ^ (1 + hihi){i + h]h,){fyi - fU]MiM - /,2/ii). 



(42) 



<i,J> 



Taking decomposition of the Heisenberg interaction term above into terms involving charge and spin fluctuations 

separately, uncorrelated mean field contributions and correlated fluctuations, i.e., correlations between charge and 
spin fluctuations as in the U(l) case, the SU(2) action is rewritten. 



S[bc 



, fa, Xi] = f dT\y^ {bldrbi^ + fldrfia) + 
•^0 i,a=l,2 



SU{2) 
t-J 



(43) 



where 



<i,3> 



+ {fi2fjl - filfj2){b]ibi2 + b\-^bj2) + c.c. 

J 



-2 E [((/i^/ji-/A/}2)(/ii/i2-/,-2/a))(l + /il/^i)(l + /i]/^,-) 

<i,j> 

+ ((1 + hlhi)il + h]hj)){fyi - fyjMlfi2 - fM 

-((1 + h\hi){i + h]hj)){{fyi - fyjMifi2 - fj2fii)) 

+ ((1 + hlhi)il + h]hj) - ((1 + h\h{){l + h]h,))) X 

((/i2/ll ^ .fll.fj2)ifjlfi2 - fj2fil) - {{fy]l - fnfj2)ifnfi2 " fj2fil)^^ 

-M^(l-/llM 

i 

- EO^^'^(/ii/i + blibi2) + iXf\fi2fa + b\^ba) 

i 

+iXf\fyii - fi2fl2 + blba - blbi2)) . 



(44) 
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VI. SU(2) MEAN-FIELD HAMILTONIAN 

The intersite holon interaction (the third term in 
Eq.(44)) is rewritten, 

((/i/li - Afh)ifnfi2 - fj2fa)){l + h\hi){l + h\hj) 



VP = -f (/i/jl - flfh){fjdi2 - fj2fil). (49) 

Here cr° is the unit matrix and cr^'^'^, the Pauli spin ma- 
trices. The fifth term in Eq.(44) is written, 



ft 



J 



J 



< |A(f > (1 + hjhi + h]hj + h]hih]hj) 



{l + h\hi){l + h\hj)) X 



(45) 



((/i2/jl - .fll.fj2)ifnfi2 - /j2/il)^ 

i(^{i + hlhi))({i + h]hj))x 



where A{j = (^.fjifi2 - fj^fii^ is the spinon singlet pair- 
ing order parameter. The intersite spinon interaction 
(the fourth term in Eq.(44)) is rewritten in terms of de- 
composed Hatree-Fock-BogoUubov channels in the same 
way as in the U(l) case. 



(/^/ji-/A/j2))((/ii/«-/i2/a 



ft ft 



2 



= ^|A/f (1+ < hlhi > + < h]hj > + < h\hi X h]hj >). 

(50) 



^{jyii - fyh)Ujifi2 - fM = vn+VE + VP, We introduced ((/^/j, - fy],){M2 - fj2fn 



ft ft 



(46) 

where vd, ve and vp are the interactions corresponding 
to the direct, exchange and pairing channels respectively. 



ft ft 



1=0 



VE 



(47) 
(48) 



(/i2/ji - flifh)){(fnfi2 - fj2fa)) = |A(f and 
((1 + hlhi){l + h]h,)) « ((1 + h\h,)){{l + h]h,)). 

To obtain the effective Hamiltonian, we introduce the 

Hubbard-Stratonovich fields, A\jafi concerned with the 
holon pairing channel in Eq.(45), and p'^, Xij and Ay 
concerned with the spinon direct, exchange and pairing 
channels. We perform the Hubbard-Stratonovich trans- 
formation for the holon pairing interaction. 



e^p("-\A{j\Hlb]0b,0b, 



cx yn„,^(iA^;„^dA^^„^exp ^-^\AQ^J2(\^iJa0 



ijaffblabif}) 



After saddle point approximation, we obtain, 

J 
2 



'■SU{2) 



(51) 



(52) 



a,/3 



with Aijajs =< biabjf) >, the saddle point of the holon pairing order parameter. 

By introducing the Hubbard-Stratonovich fields, p*^, Xij and A^ for the spinon direct, exchange and pairing order 
shown in Eqs.(47), (48), (49), we rewrite the effective Hamiltonian for Eq.(44), 



2t 



<ij> 



Q,/3 



<i,j> 



^iAifkfl - /A/I2) - T^bl,bi2 + blbj2)} - CC. 



+ tEE((^^.)'-^1(//-v.)) 



<i,j> 1=0 



9 



- 9 E l^fil' + < >-< h\hi > 



+ 



<jj> 



-E 

<i,j> 



i i 

- E (/A/i + &li«'i2) + iA?(/i2/a + blba) + iXKflifn - fi^fl + blba - blh^)], 



(53) 



where Xij 

{Mfj2 - 



^Xij;12 



/i2/,i) - j-iblb,, + bUa)) = A 



. with x-j;i2 



+ b^bii 



and x, the 



+^>^Kf}ihi-fi2fl2 + bibn 
t 

~2 



7(1^ 

hole doping rate. 

The second term in Eq.(53) represents the pairing en- 
ergy term for Cooper pair as a composite of the holon 



t ^ 
+^ E 



E 

2 



2J, 



Xij;12 



<*,i> 



pair and spinon pair. The scalar boson field, A'^-.^^ rep- ^^t ^ 

resents the holon pairing order parameter for the nearest Jp 



resents the holon pairing order parameter for the nearest 
neighbor intersite holons, namely the b^ and b^ single 
holons with the holon index, a,(3 = 1 or 2 [33]. To sim- 
plify the effective Hamiltonian, we rearrange each term 
in Eq.(53). 

For the paramagnetic state, we obtain = 
|(/./<t'/,) =< Si >= for ; = 1,2,3 and pO = i < 
fj„fi(T >= 5 for Z = 0. The one-body holon term (the 
sixth term) in the above Hamiltonian is incorporated 
into the effective chemical potential term. By setting 



6]-2^i2)(6]i&il 



bl2bj2), 



(54) 



A 



A^. 



c)Xy;l2 



and xt 



J(i 



where A{- =< fafj2 - /»2/ 



71 > 



tj:i2 — \"i2'^n ^ '^j2"^^ /' rearrange the third 
term and the second term in the bracket of the fifth term 
to obtain the effective Hamiltonian, 



<i,j> a,l3 

+^ENi^ + ^ix..r+J] 

<i,j> 

E [xUblibji - b]2bi2) - HMibi2 + b\M 

<i,j> 

-T eM.(/u>)+c.c.; 

Ai*a^(6ia6j73) + C.C. 



— CC. 



<i,j> 

- E 

2 ^ 



<i,j>.,a,l3 

-y E 

<hj> 



fj2fil) +C.C. 



- E (/A/^ + blbl,) + zXUMfii + bi2ba) 



<i,j> 
b 



where = + 2 Ej=» 

We neglect correlations between the fluctuations of or- 
der parameters. This is because correlations between 
the spin (spinon pair) and charge (holon pair) fluctua- 
tions arc expected to be small as compared to the sad- 
dle point contribution of the order parameters (the first 
and second terms in the above Hamiltonian) particularly 
near the pseudogap temperature and the bose conden- 
sation temperature. However, individual fluctuations of 
the spinon pairing and holon pairing order parameters 
are not ignored. We also neglect the fluctuations of or- 
der parameters (ninth and tenth terms) in Eq.(54). The 
ninth term represents the fluctuation of spinon pairing 
order parameter and we neglect it owing to its vanish- 
ment as an expectation value. The tenth term represents 
the fluctuations of holon hopping order parameter and is 
negligible at low temperature as its fluctuations die out in 
the low temperature regions where pairing order parame- 
ters A-^ and A** begin to open. Owing to the high energy 
cost involved with the Coulomb repulsion energy the ex- 
change interaction terms (the last positive energy terms) 
will be ignored [31]- [33]. We then obtain the mean field 
Hamiltonian, H = H^'^ + Hf + i^^ where H^'^ is the 
saddle point contribution of order parameters, Xj A-'' and 

A^ 



ttA,x 
^SU(2) 



E 

<»,i> 



Q,/3 



<»>J> 

is the holon Hamiltonian 
t 



2 , 2 
+ X 



(55) 



^SU(2) 



2 E [xiMib3i-b]2bi2) 



<i,3> 
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-A{AbUi2 + blbj2) 



— c.c. 



J 



Here the spinon quasiparticle (quasispinon) energy is 
given by 



2 J2 H\'[Al*.^p{bic.b,f,) + c.c. 

<i,j>,a,l3 

- [t^\{h\hi -x) + iX\{blb\^) + iX^ibM 

i 

and , the spinon Hamiltonian, 

J. 

<i,3> 



(56) 



with the spinon pairing energy 

= ^(jpAj^,-a(i))' + (a(2))2. 
The holon quasiparticle (quasiholon) energy is, 



(59) 



(60) 



<i,j> 



[M;{fjifi2-fi2fa)+c.c. 



El, 



(A!f)^ 



EI2 



with the holon pairing energies, 



(61) 



\f = J (A;)2 + 2\u,'>\El - 2^(M*£;^)2-(A;)2(A';-aa))2, 



= ^ (A,)2 - 2\nb\El + 2^{nbEiy-{A,nA'^-aWy 

(62) 

for the out-of-phase (the phase difference tt between 
the bi holon and 62 holon) holon pair order parameter, 



(57) 

where Xij =< flfja + f (4^ji " b]2h2) >, A^..^^ =< 
biabjfs >, A{j =< fjifi2 - fj2fii > and /if = -M - 

We take the uniform order (site-independent) param- 
eters : Xij = X for hopping, A{j = ±A/ for d-wave 
spinon pairing with the sign +(— ) for the nearest neigh- 
bor link parallel to x (y) and Aj^-.^^ = A^^ for s-wave A^^ = Ab{6a,i5i3^i ~ (5a,2<5/3,2), and the holon pairing en- 
holon pairing, with a and f3 — 1,2. For the case of ergies, 
A^^ = 0, A^*^) = and A-^ < x, the 61-bosons are 
populated at and near k = (0,0) in the momentum 
space and the 62-bosons, at and near k = (tt, tt). Pair- 
ing of two different (a ^ /3) bosons (holons) gives rise to 
the non-zero center of mass momentum. On the other 
hand, the center of mass momentum is zero only for pair- 
ing between identical(a = (3) bosons. Thus we write 
Aa/3 = Ab(5a, 1(5/3,1 ± (5a,2^/3,2) [33,35] for in-phase {+ 
sign) and out-of-phasc (— ) pairing order only between 
holons, bi and 62. Taking the uniform chemical poten- 
tial, /i'' = /i^ and the uniform Lagrangian multipliers, 
A- = A' for I = 1,2,3, the mean field Hamiltonian is 
diagonalized. 

The quasiparticle Hamiltonian is obtained to be, as 
shown in Eq.(B29) of Appendix B, 



Hl^l,^ = NJA}{2Al + X-) + NJ,(lx' + + 

k 



fe,a=l,2 



A^i = y (A;,)2 + 2\,i>'\El - 2^/(//''£;;)2 _ (a;)2(„(2)p, 
A^^ = ^J (A;)2 - 2\nb\El + 2^in»Elr - (A;)2(a(2))2 (63) 

for the in-phase (no phase difference between &i and 
62) holon pairing, A|^^ = Ab{Sa,iSp^i + (5a,2'5/3,2)- Here 

aI = JAjAb'jk, A^ = tAfipk- e{ and e\ are 
A = -tXlk, 

El = ^{ei - a(3))2 + (tAf^, - ad))^ + (a(2))2, (64) 

with 7fe = (cos/cx + cosfcy), (pk = {cos kx — cos ky) , = 
i(A(i) + A(2))/2, a(2) = (A(2) - A(i))/2 and a(3) = iX(^). 
afca(afcc) and (3ka{0ka) annihilation(crcation) op- 



+ pL^Nx. (58) erators of the quasispinons (spinon quasiparticles) and 
the quasiholons (holon quasiparticles) respectively. 



VII. SU(2) FREE ENERGY 



We rearrange Eq.(58) to obtain 
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where = NJA"^. ( 2Al 



The partition function is written, 

Z = irexp(-/3if^^) 



k,a=l,2 



+ 7V/(l + a;). 



exp(-/3if,) [](2cosh^f(l - e-''^-)-i(l - e-'^^-)-\ 



(65) 



(66) 



Prom Eq.(66), we obtain the total free energy, 

Fsui2) = NJA) {2Al + x^) + NJ' (a) + y + J) 
-2A;BT^ln(cosh(/3£;f/2)) - 2iVA;Brin2 



+kBT ln(l-e-''^-)+ 

ka=l,2 k,a=l,2 



(67) 



From the diagonaHzed Hamiltonian Eq.(58), we readily 
obtain the total free energy, 

1^ 



Fsui2) = NJA){2Al + x^) + NJ^(^A) + i 
-2A;BT^ln[cosh(/3i;^/2)] - 2NkBThi2 



X 



+kBT Y ln[l-e-^^'^'']+ Y ^^^p^+NxiiK 



fe,c«=l,2 



fc,a=l,2 



(68) 

The self-consistent equations for the amplitudes of or- 
der parameters (x, A(, and A/) are 



dF. 



SU{2) 



dx 



dF. 



SU{2) 



dAf, 



kot 

ANJAlAb 



E 

k.a 



+ 



el^Ka-l 2 



( dEf 
dx 



dEL 
dAb 



= 0, (69) 



0, (70) 



dF. 



SU{2) 



dAi 



2NJpAf + 2NAf{2A 



— j tanh 



■E 



PEf 



'dEf\ 
OA,) 

A fdEl^ 



1 2; V 



0. (71) 



The self-consistent equations for the holon chemical po- 
tential {n'') and the Lagrangian multipliers (a^'^) are 



dF} 



SU{2) 



fc,a 

, KdEl 
2^ V 



+ Nx = 0, 



dF< 



SU{2) 



- > tanh — ^ —rfr 

k 

2(e^^io. - 1) dd'') 

k.a ^ ' 



(72) 



0, fc = 1,2,3. (73) 



It can be readily proven from Eq.(73) that a^*^' = 
satisfies the three constraints above by showing that 
F(a'^^\ a^^\ a^^') is an even function of a^*^. 



f(a«)=F(0) + i^ 
+0 ((a«)^). 



2 ^ 5aWa(j) 



a(')=0 



(74) 



From the remaining four self-consistent equations in 
Eqs.(69) through (71) and (72), we determine x, A;,, A/ 
and ^ at each temperature and doping. Both T* and 
Tc arc then determined to be the temperatures at which 
the spin gap Aq and the holon pairing energy (gap) Aq 
respectively begin to open. 

As mentioned before, we consider the two cases of 
the holon pair order parameter, A^^ = A},{8a,\^0,\ ± 
'^0,2^/3,2) • There is no difference in the resulting phase 
diagram whether we choose + or — for the relative sign 
between the h\ and 62 pairing. This is because the quasi- 
particle energy dispersions become identical upto order 
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fc^ near the bottom of the quasiparticle spectrum irre- 
spective of the sign. At low temperature, the holons 
(bosons) are largely populated at and near the bottom 
of the energy. 



VIII. PREDICTED RESULTS OF PHYSICAL 
PROPERTIES FOR HIGH Tc CUPRATES 

The holon pair bosc condensation is taken care of 
by neglecting the phase fluctuations and thus by tak- 
ing phase coherence of the holon pairing order parame- 
ter. In Fig. 1 the predicted phase diagrams from the 
U(l) theory (dotted hues) and the SU(2) theory (solid 
lines) are displayed for J/t = 0.3. Encouragingly, in 
agreement with observations [8,3,10 13,18,19,21] the spin 
gap temperature (pscudogap temperature) T* tangcn- 
tially meets the superconducting transition temperature 
Tc in the overdoped region. This is attributed to the ef- 
fect of coupling between the spin and charge degrees of 
freedom. To be more specific, as shown in the pairing 
energy term, — ^\^{-\'^h\h^jhjbi for the U(l) theory and 
-i\Mj?b\j}]i^bji3h,oc for the SU(2) theory, holon pairing 
strength depends on the spinon pair (spin singlet pair) 
order which tends to disappear in the overdoped re- 
gion. The pairing energy terms here imply that bose 
condensation arises as a result of condensation of Cooper 
pairs. The Cooper pair of d-wave symmetry is a com- 
posite of a spinon pair of d-wave symmetry and a holon 
pair of s-wavc symmetry as a consequence of coupling 
between the spin (spinon) and charge (holon) degrees 
of freedom. Thus the bose condensation disappears in 
the overdoped region where the spinon pairs (spin sin- 
glet pairs) no longer exist, as observed by experiments. 

As shown in Fig. 1 the predicted pseudogap(spin gap) 
temperature Tgjj^^^ by the SU(2) theory is consistently 

SU ('2') 

higher than that of U(l) theory, Tj^^^y Tc at optimal 

doping is predicted to be lower than Tc^'^^ ■ The predicted 
SU(2) optimal doping is shifted to a larger value, show- 
ing a closer agreement with observation [8] than the U(l) 
case. Such improvement is attributed to the amplitude 
and phase fluctuations of the order parameters A-^ which 
were well taken care of in the SU(2) treatment. 

The predicted bose condensation disappears in the 
overdoping region where the spinon pair order parameter 
vanishes, in agreement with the observed phase diagram 
of high Tc cuprates. Not only the spin gap temperature 
but also the superconducting temperature is predicted to 
increase with increasing J as are shown in Figs. 2 and 3. 
This is in sharp contrast to the result of the mean field 
theories involved with the single holon bosc condensation 
[28-31,33] in which the bose condensation temperature 
(Tc txx) does not depend on the Heisenberg coupling 
energy J and increases linearly with hole concentration x. 
Thus, according to these theories the bose condensation 
temperature Tc has no correlation with the spin gap tem- 



perature T* . Indeed, the experimentally observed corre- 
lation between T* and Tc favors our proposed holon-pair 
boson theory in which the Heisenberg interaction term 
takes care of coupling between the spinon pairing order 
Af and the holon pairing order Af,. T* is determined 
from A/ and Tc from Af,. Both A/ and Ab depend on J. 
This is a cause of correlation between T* and Tc through 
the J dependence as shown in Figs. 2 and 3. However, 
a J independent universal behavior in the ratio of T* /Tc 
for the entire range of hole concentration is predicted, 
as shown in Fig. 4. The universal behavior (dotted 
{J/t = 0.2), solid {J.t = 0.3) and dashed {J/t = 0.4) 
lines) of showing a hyperbolic scaling of T* /Tc with hole 
concentration x is in qualitative agreement with obser- 
vations (shown by open and solid squares) [8,24]. This 
agreement implies that correlation between the pseudo- 
gap and superconducting transition temperatures is "me- 
diated" by the antiferromagnetic (Heisenberg) coupling 
J. 

In Fig. 5 we show the doping dependence of super- 
conducting gap at T = OJsT for both the U(l) and SU(2) 
cases. In both cases the predicted gaps decrease with in- 
creasing hole concentration. This trend is consistent with 
the ARPES [4,5] and with the Nernst effect measure- 
ments [25]. That is, a monotonic decrease of the pairing 
energy with increasing hole concentration in contrast to 
the arch (dome) shaped doping dependence of supercon- 
ducting transition temperature is predicted in agreement 
with observation. The superconducting gap is given by 
the spinon pairing energy (spin gap) Aq = 2JpAf at 
T = OK. This is because the holon quasiparticle excita- 
tion energy gap vanishes as will be shown below. 

In Fig. 6, we display the doping dependence of the 
Cooper pair order parameter, A (= A/Ab) at zero tem- 
perature for J/t = 0.2, 0.3 and 0.4 based on the U(l) 
theory. Although not shown here, the Cooper pair order 
parameter of the SU(2) theory shows a similar behavior. 
It is reminded that the Cooper pair order parameter, A 
is the composite of the spinon pairing order parameter, 
Af and the holon pairing order parameter, A^. The pre- 
dicted Cooper pair order parameter shows arch shapes 
in the plane of the Cooper pair order parameter A and 
hole concentration x as shown in the figure. This arch 
shape of the superconducting (Cooper pair) order param- 
eter was also predicted from the numerical study on the 
projected d-wave state [44]. 

In Fig. 7 we display the J dependence of the Cooper 
pair bose condensation energy U = Eisr{Ab = 0) — 
Es{Ab) as a function of hole concentration at T = OK. 
The predicted condensation energy for both the U(l) and 
SU(2) theories shows similar arch shapes as a function of 
hole concentration, exhibiting consistent increase with J 
at all doping. We omitted a plot for the U(l) case. The 
arch shaped condensation energy is in agreement with 
observation [45,46]. 

In Fig. 8, we display the ratio of the condensation 
energy to the square of the Cooper pair order param- 
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eter A = A-^A**, As shown in the figure, the ra- 

tio remains roughly eonstant for the entire range of hole 
concentration for each value of J/t. This implies that su- 
perconducting condensation energy scales linearly with 
the square of the Cooper pair order parameter, as is well 
known from the BCS theory. It is noted that this ra- 
tio increases linearly with the increase of the Heisenberg 
coupling constant J. Thus, the condensation energy also 
scales linearly with the coupling energy J, satisfying the 
relation at T = Oif , 



U ■ 



JAjAl 



JA^ 



(75) 



for both the U(l) and SU(2) theory. For the SU(2) 
theory the Cooper pair order parameter is defined as 
^Si7(2) _ y/2AfAb to take into account the contribu- 
tions of bi and &2 holons. 

Figs. 9 (a), 9 (b) and 9 (c) display the doping depen- 
dence of the superconducting transition temperature Tc, 
bose condensation energy U and superfluid weight, all of 
which show a common feature of the arch shape, again 
in agreement with observations [8,45-47]. As shown 
in Fig. 10, the predicted superfluid weight shows the 
'boomerang' behavior. This trend is in complete agree- 
ment with the well-known Uemura plot [47]. That is, the 

predicted boomerang behavior in both Tc and "° ''^^"^ 
as a function of x is consistent with the measurements 
[47,48] of muon-spin-relaxation rates a, which showed a 
linear increase of cr a ng/m* oc in the underdoped 
region and the reflex behavior for both Tc and ns/m* 
in the overdoped region. The SU(2) results are qualita- 
tively similar and are not plotted. Although not shown 
here, the predicted specific heat shows a non-Fermi liq- 
uid behavior in the underdoped region and a Fermi liquid 
behavior in the overdoped region [49]. 

It is of great interest to see whether there exists a gap 
in the quasiparticle excitation of holon (boson). The min- 
imum excitation energy for boson corresponds to the gap, 
i.e., the energy required for the quasiparticle excitation 
[41]. As will be discussed below the predicted minimum 
excitation energy at T = QK for the square lattice of in- 
finite size is found to be i?^ = O"*" at k (0, 0) for the 
U(l) theory and at k = (0,0) and (7r,7r) for the SU(2) 
theory. Excitation energy gap Eq referred here should 
not be confused with the holon pairing energy, Aq (Aq 
in Eq.(26) or Aq" in Eq.(61). We computed the holon 
quasiparticle excitation energy (gap) Eq in the thermo- 
dynamic limit N. Fig. 11 shows for iV = 20 x 20 the 
ratio of the computed Eq to the effective strength of the 
holon pairing interaction JAj at zero temperature: JAj 



comes from the pairing energy term — -^jA^j [•^6]6j6j6i for 

e'' 

for the SU(2) theory. The predicted ratio ^ ^^.^ ^ shows 

nearly a constant value of 1 /400 which is independent of 
hole concentration for the selected finite square lattice of 



the U(l) slave-boson theory and ~^\A{A'^blc^b^,fjbjpbic 



size, N 



20 X 20. As seen in Fig. 12 values of 



different values of lattice size iV, iV = 10 x 10. 16 x 16 
and 20 x 20 are independent of the lattice size iV, satisfy- 
e'' 

ing = 1. Thus the holon quasiparticle excitation 

energy at T = OK as a function of lattice size is given by 
the relation, 



JA) ^ 



(76) 



This implies that the holon quasiparticle excitation gap 
at T = QK approaches E^ = 0+ (becomes gapless (mass- 
less)) in the thermodynamic limit, N oo. Thus, 
the holon-pair bose condensed state is stable against 
the single-holon bose condensed state. This is in agree- 
ment with the theory of Nozieres and Saint James [41] 
who showed that boson-pair bose condensation is stable 
against single boson condensation since Eq > 0+ for the 
case of boson-pair [41] . 

At low temperature the holon quasiparticle excitation 
energy becomes linear in momentum. Expanding the 
holon quasiparticle energy Eq.(26) around k = (0,0) for 
the U(l) theory, one obtains the energy dispersion 



El = vk, 



(77) 



with k 



^k'^ + ky, the momentum and v = 

^2JA'jAb(2tx + JAjAb), the group velocity of the 

quasiparticle (see derivations in the previous section). 
For the SU(2) theory, the minimum of the quasiparticle 
energy occurs at k = (0, 0) and k = Q = (tt, tt). Expand- 
ing the quasiparticle energy Eq.(61) around k = (0,0) 
and k = Q respectively, one obtains the energy disper- 
sion relations, 



E', 



k,2' 



El^,=vk, 
= v\k-Q\, 



(78) 



for 



with V = ^2JAjAbitx + JAjA^) for the case of A^_^ = 
AbiSa^iSfs^i — 5a,25i3,2)- The above relation holds, also, 
true in form for the case of A^ ^ = Ab{5a,i5p,i+5a.25i3,2)- 
The linear dispersion relation with the holon group veloc- 
ity guarantees the stability of superfluidity [50]. This is 
because the quasiparticle group velocity v represents the 
critical velocity below which the superfluid flow of boson 
is stable against quasiparticle excitations [50]. In Fig. 13 
we display the doping dependence of the group velocity v. 
Both the U(l) and SU(2) theories predict an arch shaped 
group velocity in the plane of v vs. x. The predicted crit- 
ical hole doping concentration Xc at which the maximum 
velocity occurs are 0.1, 0.13 and 0.15 for J/t = 0.2, 0.3 
and 0.4 respectively. Interestingly, these critical values 
are nearly the same as the critical hole concentrations at 
which the maximum of the superconducting condensation 
energy occurs; x = 0.1, 0.12 and 0.14 for J/t = 0.2, 0.3 
and 0.4 respectively. This implies that there exists cor- 
relation between the holon (boson) group velocity and 
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the condensation energy. It is expected that the group 
velocity should decrease in the heavily overdoped region, 
since the spinon pair order parameter diminishes in 
the overdoped region. Indeed, this is predicted to be 
the case as shown in Fig. 13. As the Heisenberg cou- 
pling constant increases, the group velocity is shown to 
increase as is shown in the figure. 

The generic features of peak-dip-hump are observed in 
ARPES, STM and optical conductivity measurements, 
suggesting a common origin for its cause. Although not 
reported here, the predicted sharp quasi-particle peak 
below Tc in the spectral fimction is seen to arise as a 
result of bose condensation below Tc and its reduction 
above owing to the disappearance of the bose con- 
densation in the pseudogap (spin gap) phase [37]. This 
prediction of the sharp peak as a consequence of the 
bose condensation is in complete agreement with ARPES 
measurements [51]. Our earlier study of optical conduc- 
tivity also showed the peak-dip-hump (Drude peak-dip- 
mid-infra band) structure [38] below Tc, in agreement 
with observations [15]. We find that the hump in both 
ARPES and optical conductivity is caused by the spin 
fluctuations (of the shortest possible antiferromagnetic 
correlation length) or the spin singlet pair excitations. 
Although not shown here, the predicted temperature and 
doping dependence of spectral intensity [37] and optical 
conductivity [38] is consistent with observations. Fur- 
ther the predicted doping and temperature dependence 
of spin susceptibility is also in agreement with the INS 
(inelastic neutron scattering) measurements [14] . 

IX. SUMMARY 

In this paper, by making a full exposure to our earlier 
holon-pair boson theory [35] various physical implications 
and predicted results are discussed. We discussed the im- 
portance of coupling between the spinon (spin) pairing 
order and the holon (charge) pairing order. Supercon- 
ductivity is seen to result from the bose condensation 
of the Cooper pairs of the d-wavc symmetry which is a 
composite of the d-wave symmetry of spinon (spin) pair 
and the s-wave symmetry of holon (charge) pair. This 
theory differs from other proposed slave-boson theories 
in that coupling between the charge and spin degrees of 
freedom is manifested in the anti-ferromagnetic (Heisen- 
berg) interaction term. Further the spin-charge sepa- 
ration no longer appears in the mean field Hamiltonian 
unlike other mean field theories concerned with single- 
holon bose condensation. Both the pseudogap tempera- 
ture T* and the superconducting transition temperature 
Tc are predicted to scale with the Heisenberg exchange 



coupling strength J. Interestingly, despite such J depen- 
dence on both the T* and T^, the predicted ratio, T* /Tc 
for the entire range of hole doping concentration x is 
independent of J. Accordingly, the predicted universal 
behavior of the hyperbolic scaling T* /Tc with hole con- 
centration X is in agreement with observation [8]. The 
antiferromagnetic spin fluctuations of the shortest possi- 
ble correlation length or the spin singlet pair excitations 
are seen to induce the bose condensation; obviously the 
bose condensation does not occur in the region where the 
spin singlet pair excitations are absent as shown in the 
observed phase diagram for high Tc cuprates. The spin- 
gap temperature or the pseudogap temperature is the 
temperature where the spin-gap begins to open. Both 
the spin gap temperature and the spin gap are predicted 
to continuously decrease as the hole doping concentra- 
tion increases. Indeed, we found that the predicted bose 
condensation (superconducting transition) temperature 
drops to zero in the overdoped region where the spin 
gap or spinon pairing order disappears. This finding is 
consistent with observations. This implies that the super- 
conducting phase transition is attributed to the interplay 
between the charge and spin degrees of freedom. To put 
it otherwise, coupling between the holon (charge) pair 
and spinon (spin) pair to form the Cooper pair bosons is 
essential for causing the bose condensation. The doping 
and temperature dependence of various physical proper- 
ties are in complete agreement with observations. 

We have demonstrated that interplay between the 
charge and spin degrees of freedom are responsible for de- 
termining the characteristic features of various observed 
physical properties; the bose condensation energy, su- 
perfluid weight, specific heat, spectral functions, opti- 
cal conductivity and spin susceptibility, not to speak of 
the arch shape of superconducting transition tempera- 
ture and the tangential appearance of the pseudogap tem- 
perature. In short, using the use of the U(l) and SU(2) 
holon-pair boson theories which fully takes into account 
coupling between the charge and spin degrees of freedom 
in the antiferromagnetic (Heisenberg) interaction term, 
we demonstrated self-consistency in predicting not only 
the arch shaped bose condensation temperature but also 
various other physical properties in agreements with ob- 
servations. 



X. ACKNOWLEDGEMENT 

One(SHSS) of us acknowledges the generous supports 
of Korea Ministry of Education (HakJin Program) and 
the Institute of Basic Science Research at Pohang Uni- 
versity of Science and Technology. 



15 



APPENDIX A: DIAGONALIZATION OF THE U(l) HOLON AND SPINON HAMILTONIANS 



We first write tlie mean-field Hamiltonian in momentum space and diagonalize tlie resulting kinetic energy term. 

Wc then express the total Hamiltonian in the basis of the diagonalized kinetic energy. Finally the total Hamiltonian 
is diagonalized by using the Bogoliubov-Valatin transformations. We take the uniform, that is, site-independent 
hopping order parameter [31], Xij = X^^'^, where the sign +(— ) is for the counterclockwise(clockwise) direction 
around a plaquette and the uniform pairing order parameters, A^- = A/e^*"^ and A^^ = Afce^*'^ , where the sign 
+(— ) is for the ij link parallel to x (y) where A^, A/ and x a-re the site-independent amplitudes corresponding 
to proper order parameters above. Wc also take the site-independent chemical potentials, iJ,j — and /xf = jJ' 
[29,31,34]. 

The saddle point energy term in Eq.(20) is, then. 



J 



A 



<i,j> 

--nj[a}aI + a 



+ NJ, 



^ +4. 



(Al) 



with A'', the lattice size. 

The kinetic energy term of the holon in Eq.(22) in momentum space is obtained to be 



ix*jtb]h + c.c.) 



<i,3> 



N 

k,k' <».J> 



Xlk cos 6 
-iX'Pk sin( 



iXfk smfc) 
-X7feC0s( 



h+Q 



(A2) 



and the kinetic energy term of the spinon in Eq.(23) in momentum space, 



H 



K.E 



4 N 



k.k' <i,3>,(' 



9 JZ'- -l^ka- fk+Q,a ) 



k(T 



Xlk cos 6 
-iXfk sin( 



iX(pk sm 6 
-Xlk cos ( 



fka 
fk-\-Qa 



(A3) 



where 7^ = {coskx + cos ky) and (pk = (cosfc^ — cosfcj,). 
The holon pairing term in Eq.(22) is 



K = -J^ ^[a5:(&a)+c.c.] 



<i,J> 



= -J- 



2 N 



k,k' J 



-J-^Afc^[(cosT^7fe + i sin T''(pk){bkb-k) + ex.], 



(A4) 



and the spinon pairing term in Eq.(23) in momentum space. 



J, 



16 



2 



k.k' 3 



= — ^A/^[(cosT-^7fe +isinT-^(^fc)(c7/fe<^/_fe_cr) + C.C.]. 

k 

The chemical potential terms of spinon and holon respectively are 

VE(^l^^-^) = VE^I^'»+'"'^^- 



Finally, combining Eqs.(A3), (A5) and (A7), we write the spinon Hamiltonian, 



k<7 



X7fc COS ixfk sm & 
-iXifkSine -xik cos 



fka 
fk+Qa- 



J, 



f^f^{i^°^'^^'yk+is^^'^^V'k){(^fkaf-k-a) +C.c: 



VE/L/fe-+/^(i-^)> 

ka 



and combining Eqs.(A2), (A4) and (A7), we write the holon Hamiltonian 



Xlk cos ixfk sin 9 
-iXipk sine -Xlk cos e 

k 

-J-^Ab ^[(COS tSa; + «sin TVfc)(^'fc^-fc) + C.C.] 



hk 
bk+Q 



-li''J2bibk + fJ-^Nx. 



(A5) 

(A6) 
(A7) 



(A8) 



(A9) 



From Eq.(A3) and Eq.(A2), the kinetic energy terms 
of the spinon and holon parts are rewritten, 



^K.E. — o 



k,a ^ 



(AlO) 



(All) 



with 



Nk{x,0) 



Xlk cos 6 ix^Pk sm ti 
-iXipkSin0 -xikCosO 



In order to diagonalize the kinetic energy term, we 
introduce a basis transformation between bare particle 
(holon 6 and spinon /) basis and quasi-particle {B and 

F) basis, 



( Bk+\ _ 



= ul 



bk 



(A12) 



with U 
with 



Fka+ 
Fka- 

Uk+ Uk- 
Vk+ Vk- 



a unitary transformation matrix 



Uk±\ _ 1 f xikCose±x£.k{d) 



Vk± 



N± 



-iX'fk sm( 



(AM) 



and 



N± = V2xaW(xaW±X7fcCos^). (A15) 



In the quasi-particle basis, the kinetic energy term is 
rewritten. 



hLe= E ^xsik{e)Fl,Fu. 

/c,cr,s=±l 

H'k.E= E '^tXsUQ)BlBus, 

/c,s=±l 



(A16) 
(A17) 
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with ^k{0) = cos2 6* + V5| sin'^e. Here, s = 1 (-1) for holon. 

represents the upper (lower) band of the quasiparticle en- Fi'^ally, inserting Eqs.(A16) and (A18) into Eq.(A8) 

ergy spectrum in the reduced Brillouin zone defined by ^nd Eqs.(A17) and (A19) into Eq.(A9), we obtain, re- 

I'i'xl + l^i/l < and J2 > the summation over the reduced spectively. 



Brillouin zone 



The chemical potential terms in the quasiparticle basis „/ _ / / /M?t i? _i_ / ^T^^ \ 

are, ^K.E.+chem. — \^ks ~ )^kas^kas + fJ- {i- - X) 

k,(T,S—±l 

fkafka- = -^-^ {jkafka + fk+Qafk+Qaj ^ ^ 

, tor spmon and 

= -M^^(f,V+F,.+ +F,V_F,._) (A18) 

, . ^ Hl^.+chem.= ^'ks-l^')BlBu, + ,J>Nx (A21) 
tor spmon and 

— /Lt^ b\hk = —fj!' {^\bk + b^+Q^k+Q^ for holon. Here the spinon quasiparticle ( "quasispinon" ) 

k k energy is = ^xsi,k{d) and the holon quasiparticle 

hv-/r.t r.\ r. \ . A , ^^ ( " quasiholon" ) energy , e^^ = 2txs^k{6) with ^ (6') = 

I The spinon pairing term in Eq.(A5) is rear- 



ranged to yield. 



Bp = --yA/ ^[(cosT-''7fe + ismT^ ipk){afkaf-k-a) + c.c] 



2 

k 



= -YA/^[(cosT-''7fe +isinT-^(pfe)(^(/fe^/_fe^ + f-k^fki) 

k 

-{fk+Q^f-k-Ql + f-k-Q^fk+Qi)) + C.C.]. (A22) 

Using the symmetry 7_fc = 7fe and v'-fe = <Pfe, we rewrite, 

-f^p = - Jp^f X^[(cosr-''7fe + is\nT^ip},) {fk\f-ki - fk+Q\f-k-Qi) + c.c] 

k 

= -Jj,AfY}ak{Tf){Pk^+F-ki+ + Fkf_F_k^_) + c.c], (A23) 



where ak{T^) = {cosr^'^k + ismr^ ipk) = e'*^fc(T^) with ^k{T^) = ^ 7^ cos^ + ip^ sin^ and tan$ 

that is, 



^ tanr-''. Here we used the unitarity of 



Uk+ Uk- 
Vk+ Vk- 



UksU-ks - VksV-ks = {uksT - {vksY = \uks\^ + \vks? = 1 

UksU-k-s - VksV^k-s = {Uks)"^ - {Vk-sf' = UksU*k_^ + VksV*k_^ = 0, (A24) 

and u*^^ = Uks and v^^ = —Vks for s = ± from Eq.(A14). 

Similarly, the holon pairing term in Eq.(A4) is rearranged to yield, 

Hi = -^lA/I^Afc VKcosAfe +«sinrVfc)(^fc^'-fe) +C-C.] 
^ k 

= -^\Af\^Ab Y,[ak{T''){Bk+B-k+ + Bk-B-k-) + c.c.]. (A25) 

k 
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Prom Eqs.(A20) and (A23), we obtain the effective Hamiltonian in the transformed basis, 

= E (^L - f^^)FLPk.s - ^[e^''^^^">AlF,^sF-kis - c.c] + /.^iV(l - x), (A26) 

k,s,<T k,s 

for spinon and from Eqs.(A21) and (A25), 

H''= i^l-t^'MsBke- E [e'''''^^'^^BksB_ks + c.c.]+f,''Nx (A27) 

fe,s=±l k,s=±l 

for holon. Here = ^sx^kiO), U^) = cos'' + sin'' 9 , = JpAf^kir^), = 2tsxUe) and = 

Taking Fkas = e-'^'-^'^^^/'^F'^^^ for spinon and Bks = e-**''^^'')/^^^^ for holon, we rewrite Eqs.(A26) and (A27) 
respectively 



= E ('L - f^^)F'kisF'kas - E ^liKs^F-ksi - f'^u^-Hs) + - x), (A28) 

fc,s,cr k^s 
I I 

H''= E E %B',,B'_,,+BlB\^)+^,''Nx. (A29) 



fc,s=±l k,s=±l 



For simplicity, we will omit the prime symbol ' in Ff.^^ and B^.^. 

Introducing the Bogoliubov-Valatin transformation for spinon quasiparticle operators aksa- and cx].g^, 

Fkts = cosCfesQ!fesT +sinCfesa^^fesp 
F-kis = - sin Cksai^-^ + cos (ks^-ksi , (A30) 

we rewrite Eq.(A28), 

k,s k,s 



k.s 



(A31) 

where is the spinon quasiparticle (quasi-spinon) energy, E^^ = \J (ejj^ — ji^Y + (Aq)^ and Aq is the spinon pairing 

energy (gap), A^ = JpAf^kir^)- 

Introducing the Bogoliubov-Valatin transformation concerned with holon quasiparticle operators pks and f^l^, 

Bks = coshr;fe5/3fes + sinh?]^^/?!^,^, 
B-ks = smhriksPls + coshijksP-ks, (A32) 

we rewrite the holon Hamiltonian, 

H'= Yl Fljlf3ks+ E l{El+f^'')+f^'Nx, (A33) 



fe,s=±l fe,s=±l 



where E^.^ is the holon quasiparticle (quasi-holon) energy, E^g = (e^j, — jjf'Y ~ i^oY ^'^d Ag is the holon pairing 
energy, Ag = JAjAt^k{T^). In the last term we used the identity J2k,s=±li~^ks+^^'') = ~ Sfe,s=±i(2*sxCfe(^)-M'') = 



Ek.sf^'- 
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Combining Eqs.(Al), (A31) and (A33), we obtain the total quasiparticle Hamiltonian, 

1 o 1- 



A2 + -v2 ^ - 

f ^2^ 4 



I I 



fe,s=±l 



fe,s=±l 



(A34) 



APPENDIX B: DIAGONALIZATION OF THE SU(2) HAMILTONIAN 

Here we diagonalize the SU(2) Hamiltonians for both the spinon and holon in Eqs.(57) and (56). For the saddle 
point energy (the second term) in Eq.(55), we write 



Si7(2) 



<i,j> a,/3 <i,j> 



= NJA}{2Al + x^) + NJp{A} + + \), 
where we used A'-j.^^ = Ab(^a,i^/3,i ± ^0,2(^/3, 2)- The spinon Hamiltonian in Eq.(57) is rewritten, 



(Bl) 



(B2) 



-2A/, 



where Vi = ( j and Uij = J J^^^^ ) with a^^) = a^^) = -^'+^' , a^^) = iX^ and f, the Pauh 



spm matrices. 

For the uniform phase of the hopping order parameter and the d-wavc phase of the spinon pairing order parameter 
( Xij = X) M,i+x = ~M,i+y = A/ ), the order parameter matrix Uij is given by, 



i,i+y 



X -2A/ 
-2A/ -X 

X 2A/ 
2A/ -X 



(B3) 



Using Eq.(B3), we rewrite the spinon Hamiltonian in momentum space. 



H: 



f 

SU(2) 



J, 



<i,j> 



JpA/^,-aW- 



(2) 



(2) 



-fXlk + a 



(3) 



(B4) 



where tpk = { ] and 7^ = (cos kx + cos ky) and = (cos — cos /cy). 

/-/s2 



The spinon Hamiltonian in Eq.(B4) is rewritten. 



where e{ = -^xik and A{ = JpAf^pk - a^^) To 
diagonalize the Hamiltonian, the unitary (Bogoliubov- 
Valatin) transformation is introduced. 



(B5) 
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fkl 
k2 



/I 



aki 
t 

k2 



a 



(B6) 



H: 



SU(2) 



-Ei 



and ?;/ = 







with ^ ^. In the quasiparticle basis, the spinon with eI = ./I4^^^^¥h^^ = -J.X.k and 



Hamiltonian (B5) is diagonaUzed, 



= |A{| = V(JpA/<pfe-a(i))2 + (a(2))2. 

The holon Hamiltonian in Eq.(56) is rewritten, 



-"5(7(2) 



4 I] [h\Uijhj + c.c] 



<hj> 



J 



-^/iJ(ai-r + M'')/ii 



+Ij!'Nx, 



(B8) 



where /ij 



■ 



and f7,j- = 



/* 



with a(i) = a(2) = ~^^+^^ , a^^) = iAf and f, the Pauh spin 



matrices. 

Taking the uniform phase of the hopping order parameter [xij = x) and the d-wave symmetry of the spinon pairing 
order parameter (AJ_._,_^ = —A{-_^_y = A/), we write, 



TT. -fx -A/ 



X A/ 
A/ -X 



(B9) 



Using the above expression in Eq.(B9), we rewrite the first and third terms of Eq.(B8) in momentum space, 



<i,3> 



(3) _ 



)hk, 



(BIO) 



where /ife = ^ ^'^^ J and 7/5 = (cosfca, + cosky) and t^fe = (cosfcx — cosfcy). 

For the s-wave symmetry of the holon-pair order parameter (A^_,_^ = A^_^_y = A^) with the same ampUtude for 
the 61 — bi and 62 — ^2 holon pairs, the holon-pair order parameter matrix is written. 



A^ = 



Ab 

e"f'Ab 



(Ell) 



where e"^ is the phase difference between the bi — bi and 62 — &2 order parameter. Here we consider the two cases of 
_ _|_]^_ Transforming the holon pairing term (the second term of Eq.(B8)) into momentum space, we obtain 

J 



^P = -I E \^{,\'[hlAfj{hlf + c.c. 

<i,j> 

= -^iA/PEh4A^('^-.r+- 



(B12) 



with 7fe = (cosfca; + cosky). Here we used the identity h_f.A^{h 
hl{A^)^{hl^)^ = hlA^ihl^f and (A^)^ = A^, where (A^)^ is the transpose of A^ 



°-ka^a(i'^k0 



b^ b'^ 
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Inserting Eq.(BlO) and Eq.(B12) into Eq.(B8), we obtain the holon Hamiltonian in the momentum space 

tisU{2) - 2^H\ f^^^^ _ ^(1) _ ,;„(2) ^ „(3) _ „fe ) 



tXlk + 



-^iA,i^^[7.4A^(/.L,r+c.c; 

k 

+IJ^'Nx. 

The above holon Hamiltonian is rewritten, 

k 

= E \hlMkhk + h}_^M_kh-k 

k 

+ E [hlNkihlJ^ + hl^Nlhk] + fJ'Nx, 



+ n''Nx 



(B13) 



(B14) 



where 



,(3) 



2 V tAfipk - a(^) - m(2) 



,(2) 



)■ 



(B15) 



In order to represent the above holon Hamiltonian in terms of the 4-component holon field operator, we rearrange 
the second term of Eq.(B14), 



M .h (b^ )( ^->^ ] ( ^-^1 



= M]_\bl^,b^ki + 6Lfei6-.2 + Mf_ibl,,b.ki + Mf_ib\^^b. 
Inserting Eq.(B16) into Eq.(B14), the holon Hamiltonian leads to 

k 

+ E [hlmhlkf + hlkNlhk] + n'Nil + x) 

k 

WHlMlHk+n''N{l+x), 
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f22 .t 



-k2 



-"5(7(2) 



(B16) 



(B17) 



where Hk 



( bki ^ 

bk2 

'-kl 

\ blk2 J 



and 



tA 



ad) 



-tXlk - fJ." -a 



(3) tAftpk- a'^^^ +ia 



7,(2) 



txik - 1^" + a 



-JAjAll^k 
-JA2A5*27k 



(2) -JA2A5i7k -JA2A^27k 

(3) -JA2A^i7k -JA2A^27k 

JA2A|*i7k -ix7k - - a(') tA/c^k - a« - 

JA2A||7k iA/(^k-ad)+«a^^) tX7k - M*" + a^^^ 



(B18) 



with 7k = (cosfcx + cos ky) and ipu = {coskx — cos ky). 
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To diagonalize the above Hamiltonian, we consider the 
Bogoliubov-Valatin transformation for the hofon field op- 
erator iJk, 



where Ea is the quasiparticle energy. The Eq.(B21) can 
be rewritten, 



(B19) 



[Ml - E^J] • wg = 0, 



(B22) 



where H^(k) is a 4 x 4 matrix and H^, hofon quasiparti- 
cle ( "quasiholon" ) operator. For the quasiholon operator 

to satisfy the boson commutation relation, the trans- 
formation matrix VF(k) should satisfy 



where J 



W^(k)JW^(k)^ = J, (B20) 
It is noted that VF(k) is 



10 
10 

0-10 
-1 

not a unitary matrix unlike the Bogoliubov-Valatin trans- 
formation for fcrmion operators. VK(k) is chosen so that 
the transformed Hamiltonian matrix H^(k)^^M^H^(k) is 
diagonalized, that is, 



W"(k)tM^W^(k) 



(El 











\ 





E2 














—E3 







V 








-E4 


J 



(B21) 



where is the a-th column vector of W{k) with a = 
1, 2, 3, 4. The quasiholon energies and the corresponding 

eigenvectors are obtained from the above equation. This 
equation is solved by the symbolic calculation of Mathe- 
matica [52]. 

Only pairing between the same kind of holons allows 
the pairing state of the zero center of mass momentum. 
At low temperature bi and 62 holons are equally pop- 
ulated at k = (0,0) and (7r,7r) respectively. Thus, we 
consider the same amplitude for boson pairing and 
b2-b2 boson pairing, but not for 61-62 boson pairing. To 
allow the phase difference between the pairing order pa- 
rameters of these two channels, we consider the two cases 
of A^^ = Ab{Sa,iSp,i ± (Ja, 2(5/3, 2) with + sign for no 
phase difference and — sign for the phase difference tt 
between the bi and 62 boson pairing. For the case of 
AJ^ \ / 1 

eigenvalues are given by 



the quasiparticle energy 



El = y {E^,r + i^^'Y - W + 2^1^)^ - (a;)2(a; - J^, 

EI2 = sj {Elf + - (a;)2 - 2.J{i^Y - (a;)2(a; - 

El, = -^{Eir + iii^r - (a;)2 + 2^/{iJ^if - (a;)2(a; - 5^ 



El, = -J {Elf + {f,ff - (a;)2 - 2^i„^Elf-{A',f{A';-a(^)f 



(B23) 



A?2A 

A^2; 








El = 


^iElf + {f,»f- 


{A',f + 2^{^,'>Elf- 


{\f{a(^)f, 


EI2 = 


^{Elf + {n^f- 


{\f-2^{l,^Elf- 


(A;)2(a(2))2, 


Els = 


-^{Elf + in^f 


-{A',f + 2^i,j,^Elf 


- (A;)2(a(2))2 


Eli = 


-^{Elf + if.'^f 


-{A',f-2^{fj,bElf 


- (A,)2(a(2))2 



(B24) 



where El = ^(-iX7k)2 + (A^(k))2, A^ = tAfcpu, and A^ = JA2Afe7k. 



23 



Inserting Eq.(B23) or Eq.(B24) into Eq.(B21), we ob- 
tain 



Vr(k)tM^W(k) = 



( El^ Q \ 
El^ Q 

Q El^ Q 



(B25) 



V Q El^J 
and the diagonalized holon Hamiltonian, 

Hsui2) = \Y. HlMtH^ + n'Nil + x) 

k 

= J E H'^W{\,)^M^W{\<.)Hi, + t,'N{l + x) 
^ k 

^ k 

+Ei^{pyk2 + Pk2f3i^)] + l^'Nil + x), (B26) 

where we used the fact that E^^ is the even function 
of k. After normal ordering, we obtain the diagonaUzed 
Hamiltonian, 

Hsui2) = E [^LPijka + + + x), (B27) 



k,a=l,2 



where 



(B28) 



Combining Eqs.(Bl), (B7) and (B27) we obtain the 
diagonalized total Hamiltonian, 

H^Jl,^ = NJA}{2Al + x') + NJ,[^x' +^} + \) 
+ E-^fc("fei"*i ~ "fe2aj[.2) 

k 



fe,a=l,2 



fi'Nx. 



(B29) 
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^ 0.07 




X 

FIG. 1. Computed phase diagrams with J/t = 0.3. 
Tgjj^2) {Tu{i)) denotes the pseudogap temperature and T^^^'^^ 
(Te^'^^), the holon pair bose condensation temperature pre- 
dicted from the SU(2) (solid lines) and (U(l)) (dotted Unes) 
slave-boson theories respectively. The scale of temperature is 
based on t = 0.44ey[43]. 
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FIG. 2. Doping dependence of the holon pair bose con- 
densation (superconducting transition) temperature Tc for 
J/t = 0.2, 0.3 and 0.4 based on the U(l) and SU(2) holon-pair 
boson theories. 
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X x/Xp 

FIG. 3. Doping dependence of the spin gap (pseudogap) FIG. 4. The ratio of the pseudogap temperature to the 

temperature T* for J/t = 0.2, 0.3 and 0.4 based on the U(l) superconducting temperature, T* /Tc as a function of scaled 
and SU(2) holon-pair boson theories. hole concentration. Experimental ratios (squares) arc for 

Bi2Sr2CaCu20s+6 and La2-xSrxCu04[S,21]. Xo is the opti- 
mal doping rate. 
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FIG. 5. The doping dependence of the superconducting gap 
= 2JpA/ at T = for J/t = 0.2, 0.3 and 0.4 based on the 
U(l) and SU(2) theories. 
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FIG. 6. Doping dependence of the Cooper pair order pa- 
rameter, (A/Ab) at T = for J/t = 0.2, 0.3 and 0.4 based 
on the U(l) theory. 
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FIG. 7. Doping dependence of the condensation energy U 
at T = for J/t = 0.2, 0.3 and 0.4 based on the SU(2) theory. 
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FIG. 8. The ratio of the condensation energy, U to the 
square of the Cooper pair order parameter A'^ as a function 
of hole concentration at T = for J/t = 0.2, 0.3 and 0.4 
based on the SU(2) theory. 
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FIG. 10. Superfluid weight, "^C^^^^ooit) superconduct- 
ing transition temperature, Tc in the U(l) theory with 
J/t = 0.2. The directions of thick arrows denote increasing 
hole doping concentration. Xo = 0.07 is the predicted optimal 
doping rate. 
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FIG. 9. Doping dependence of (a) superconducting tem- 
perature, (b) the condensation energy and (c) the super- 
fluid weight with J/t = 0.2 for the U(l) theory. The su- 
perfluid weight and the condensation energy was obtained at 
T/t = 0.001 (equivalent to 5K with the use oit = 0.44ey [43]). 
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FIG. 11. Doping dependence of -j-^, the ratio of the holon 

quasiparticle excitation energy gap Eq to the effective holon 
pairing interaction strength JA^ at T = for J/t = 0.2, 
0.3 and 0.4 using the U(l) and SU(2) theories respectively. 
AT = 20 X 20 lattice is used. 
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FIG. 12. Doping dependence of N—-^, the holon quasipar- 

f 

tide excitation energy gap Eq to the effective holon pairing 
interaction strength JA^ multiplied by the lattice size A'^ at 
T — using the U(l) and SU(2) theories respectively. Three 
different lattice sizes of A'^ = 10 x 10, 16 x 16 and 20 x 20 are 
used with the choice of J/t = 0.2. 
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FIG. 13. Doping dependence of quasiholon velocity at 
T = based on the U(l) and SU(2) theories respectively. 
AT = 20 X 20 is used. 
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